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1. Introduction 

The AdS/CFT correspondence provides a precise map between observables of a quantum 
field theory at strong coupling, and classical fields in a weakly coupled gravitational theory. 
Since the gravitational theory lives in a spacetime with one dimension more than the field 
theory, the correspondence is also known as a holographic duality. This map provides 
an obvious advantage for computing quantities at strong coupling, and has motivated the 
study of many toy models aimed at applying the duality to real systems - traditionally to 
QCD but also more recently to strongly correlated condensed matter systems (see [1-4] for 
reviews on these topics). 

One topic of interest for condensed matter applications of holography, or AdS/CMT, 
is the realization of Josephson junctions [5]. These objects consist of two superconductors 
separated by a "weak link" made of another material: an insulator for "SiS" junctions, a 
normal metal for "SnS" junctions, or simply a narrowing of the contact surface for "SsS" 
junctions. Generically, Josephson junctions describe interfaces in superconductors across 
which the electron-pair condensate suffers a change in phase. For SiS and SnS junctions, 
quantum mechanical tunnelling across the interface induces a non-zero current proportional 
to the sine of the phase difference, even in the absence of applied voltage. This is the DC 
Josephson effect. Meanwhile, if one applies a non-zero DC voltage across the junction, one 
observes a current which oscillates in time (the AC Josephson effect). 

A number of works have recently used AdS/CMT techniques to model Josephson 
junctions in various dimensions and configurations, for s-wave [6-9, 11] and p-wave [10] 
superfluids, including a grid of Josephson junctions in [12]. These works are based pri- 
marily on the Abelian-Higgs model for holographic superfluids, with boundary conditions 
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that break translational invariance. While quite fruitful for numerical studies of the phase 
structure, these models typically require one to solve complicated coupled partial differ- 
ential equations. The aim of this work is to construct a simpler geometrical picture in 
terms of a D-brane model. In particular we will present what to our knowledge is the first 
explicit realization of a non-Abelian junction in holographic models. In principle a non- 
Abelian Josephson effect could appear in systems with an order parameter charged under 
a non-Abelian global group, such as superfluid s He [14], high T c superconductors [15], 
Bose-Einstein condensates of atoms with spin [16] or the CFL (Color-Flavor-Locked) su- 
perconducting phase of QCD [17]. 

The theories we study are (p + l)-dimensional super symmetric U(N C ) gauge theories 
at strong coupling, with a large value of N c , and correspond to the low-energy theory living 
on a stack of Dp-branes. These theories have an SO(9 — p) global symmetry group, that 
is partially broken when a few branes are separated from the stack. 1 In this sense the 
system is dual to a superfluid phase. The magnitude of the symmetry-breaking condensate 
corresponds to the radial distance by which a small number of Dp-branes are separated 
from the 0{N C ) stack. 

To this holographic p-dimensional superfluid, one can add a co-dimension one defect 
by including a D(p + 2)-brane intersecting the Dp-branes. The defect explicitly breaks the 
global SO(9 — p) symmetry to a SO (3) x SO(6 — p) subgroup, which is then spontaneously 
broken to SO (2) x SO(5 — p). When the N c — > oo limit is taken, the theory on the 
Dp-branes has a holographic dual description, with probe Dp and D(p + 2)-branes in a 
background geometry. These theories were first considered in a holographic context in [18]. 

The defect is such that the (non-Abelian) phase of the condensate can be different 
on either side of it. Therefore, the Dp/D(p + 2) intersection can be interpreted as a 
holographic realization of a Josephson junction. The effective description of the brane 
construction is very similar to the field theory model of a non-Abelian junction proposed 
in [13]: the fluctuations of the Dp branes on both sides of the defect are the Goldstone 
bosons of the spontaneously broken symmetry. In this way the symmetry is naively doubled, 
but interactions between the fields on either side break the symmetry to the diagonal 
subgroup, so one is left with the correct number of true Goldstone bosons, while the other 
modes become pseudo-Goldstones and are responsible for the Josephson effect. In our case 
supersymmetry prevents the appearance of a DC Josephson effect as we discuss in more 
detail below, but additional fluxes can be turned on in the D(p + 2) brane, that induce an 
AC Josephson effect. 

Though the Dp/D(p + 2) brane intersection is certainly an idealized model, it never- 
theless shares some interesting properties with the physical examples mentioned above. In 
particular the ground state is degenerate, meaning that there is a global symmetry that 
remains unbroken in the superfluid phase. Let us assume that there are some parameters 
that determine the properties of the junction, and that these parameters can be changed 
adiabatically, in such a way that the evolution forms a closed curve in parameter space. 
Since the ground state is degenerate, for a non-Abelian junction the final state of the 

lr The gauge group is also partially Higgsed. 
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system does not necessarily coincide with the initial state: they may differ by a Berry 
holonomy [19]. 

In the Dp/D{p + 2) intersection the parameters of the junction are magnetic fluxes on 
the D{p + 2) brane, and the Berry holonomy depends on the amount of electric flux, or 
equivalently, on the number of Fl strings dissolved on the D{p + 2) brane. One can see the 
Berry holonomy thus defined as a topological property of the intersection, and it can be 
measured through the Josephson current on the Dp brane induced by the magnetic fluxes. 
Other definitions of Berry holonomies in brane intersections are also possible: for example, 
they have been studied in black hole systems corresponding to D1/D5 intersections [20] 
and for pairs of DO branes circulating around each other [21]. 

The D3/D5 intersection deserves special attention as it was studied in detail by Gaiotto 
and Witten in the context of supersymmetric boundary conditions and S'-duality [22- 
24], including fluxes on the D5 brane [22]. Even with the string flux on the D5, the 
system remains supersymmetric. The effect of the flux can be seen as a modification of the 
boundary conditions for the fields living on the D3 brane. Although this goes beyond the 
scope of this paper, it would be interesting to see the relation to the Berry holonomy. 

The paper is organized as follows: in Section || we present the Dp/ D(p+2) intersections 
with electric flux and discuss their holographic interpretation as models of supersymmetric 
superfluids with a defect. In Section ||] we show how magnetic fluxes on the D(p + 2) brane 
induce a Josephson current on the Dp worldvolume in the Abelian case. We generalize to 
the non- Abelian case in Section |], and we compute the Berry holonomy and connection in 
Section [|. We end with the conclusions and future directions of the results in Section |(| 

2. Supersymmetric superfluids with a defect as Josephson junctions 

In this section we will present some simple models that describe supersymmetric superfluids 
in different dimensions. The superfluids have a co-dimension one defect where a Josephson 
effect can be induced by changing some parameters. The effect is non-Abelian, and it can 
be characterized by a Berry holonomy that we will define in section ||. 

Our starting point is the following D-brane configuration, consisting of a D3 and a D5 
with N s units of electric flux E on the D5 worldvolume: 



# 




0123456789 


N c 


D3 


• • • • 


1 


D5 


• • • • • • 


N s 


E 


• • 



(2.1) 



There are iV c D3 branes intersecting a D5 brane 2 localized in the 3 direction. The D5 with 
electric flux can be interpreted as a bound state of a D5 brane with N s Fl strings. This 
is a 1/4 BPS configuration, as one can see by taking T-dualities along the 1256 directions, 
which results on a configuration consisting of an (N s , 1) string ending on a D3 brane. 

One can easily generalize this construction for generic Dp/D{p + 2) intersections. The 
low-energy theory on the D3 branes is 3 + 1 dimensional N = 4 U(N C ) super- Yang Mills. 

2 The setup can be generalized further by incorporating Nf coinciding D5 branes. 
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T-dualities along the 1 and 2 directions also allow one to construct lower dimensional 
supersymmetric theories: a 2 + 1 theory on D2 branes intersecting a DA brane and a 1 + 1 
dimensional theory on Dl branes intersecting a D3 brane. One can also generate higher 
dimensional theories by taking T-dualities on the 7, 8 and 9 directions, giving D4/D6, 
D5/D7 and D6/D8 intersections. 

The Dp branes can end on the D{p + 2) brane, so the Dp branes can be split in two 
halves and separated along the directions parallel to the D(p + 2) brane. Prom the point 
of view of the fields of the D{p + 2) brane, the Dp/D(p + 2) intersections carry magnetic 
charges. From the point of view of the Dp, the intersection with the D(p + 2) brane is a 
codimension one defect. 

We can construct a holographic dual of a superfluid phase for p < 4 as follows. We 
start taking the number of Dp branes N c to be very large. In this limit, one can substitute 
the Dp branes by the geometry they source, and the Dp low-energy theory is described 
holographically by the near-horizon geometry [28]: 

^ = (f) (7 "' ,/2 ^ + (^) M/2 ^ W + tf(^) (p - 3,/2 <i nt p . (2.2) 

Here dUg_ p is the metric of a unit (8 — p)-sphere, S 8 ~ p . The "conformal boundary" is at 
r — > oo. Except for p = 3, there is a non-trivial dilaton profile 

/ p\ (7-p)(3-p)/4 

e<f = efo f^j . (2.3) 

The only other background field that has non-trivial profile is an electric (p + 2)-form flux 
on the directions transverse to the sphere. 3 The potential F p+ 2 = dC p +± is 

C p+1 = — (L) 7 ~ P dr A dx° A • • • A dx p . (2.4) 
9s 

The flux equals the rank of the group of the dual field theory or equivalently the number 
of Dp branes. The D{p + 2) brane becomes a probe in this geometry, wrapping a two 
sphere in the S s ~ p , and extended along time, the radial direction and p — 1 of the spatial 
directions x l . 

For p = 3, the geometry is dual to 3 + 1 dimensional Af = 4 SU (N c ) super Yang- 
Mills. One can obtain the lower dimensional theories through dimensional reduction. The 
field content includes six real scalars a = 1 — 6 in the adjoint representation whose 
eigenvalues parametrize the Coulomb moduli space that is locally M 6 . For instance, if one 
of the eigenvalues of one of the scalars has an expectation value, 





V ■■•/ 



(2.5) 



3 In the p = 3 case the F& flux should be self-dual, so there is also flux on the S* 5 . 
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the gauge group will be Higgsed U(N) — > U(N— 1) x U(l) and there will be an spontaneous 
breaking of R-symmetry 50(6) — >• 50 (5). 4 Therefore, states on the Coulomb branch of 
N = 4 can also be seen as (non-Abelian) superfluid phases, with 6 x N possible order 
parameters (scalar condensates) that describe the breaking of both global and local sym- 
metries. The story is similar for the lower dimensional theories, but the number of scalars 
increases by the difference in the number of dimensions, so the R-symmetry group is 50(7) 
in 2 + 1 dimension and 50(8) in 1 + 1 dimensions. 

In the brane picture, giving an expectation value to the scalars amounts to separating 
some Dp branes from the stack. If this number is very small compared with the total 
number of branes, the holographic dual description involves introducing a small number of 
Dp branes probes in the dual geometry and neglecting their backreaction. These branes 
are localized in the radial and the sphere directions, and extend along the coordinates. 
The position in the radial direction is proportional to the absolute value of the condensate, 
while the position on the sphere fixes the breaking of global symmetry, as the original 
50(9 — p) isometry of the sphere is reduced to a smaller subgroup. 

We call an idealized Josephson junction in three spatial dimensions a two-dimensional 
surface that connects two regions of space in which the condensate has a different phase. In 
other words, the spontaneous symmetry breaking is different on either side of the junction, 
even though the absolute value of the condensate may be the same. In our language, it 
should correspond to a defect where a Dp brane changes its position in the sphere as it 
crosses it. Such configurations can be realized by adding additional D(p + 2) branes in the 
way we have described before. The intersection between the Dp branes and the D(p + 2) 
is (p — 1) + 1 dimensional and will correspond to the surface of the Josephson junction. 
The defect dual to the D(p + 2) brane carries its own degrees of freedom. The field theory 
on the defect has a global U(l) symmetry, 5 for which one can introduce a nonzero charge 
density. 6 

The Dp brane can end on the D{p + 2), so instead of having a single Dp brane in- 
tersecting the D{p + 2) brane, one can "break" the Dp brane in two halves and move the 
endpoints of each half to different positions on the D(p + 2). In this case, the endpoints can 
move in the radial direction and in the directions of the 5 2 that the D(p + 2) wraps. In the 
first case the absolute value of the condensate on each side of the defect is different, while 
in the second case the pattern of symmetry breaking (the phase) is different. In particular, 
if we move the Dp branes along a maximal circle in the 5 2 , leaving everything else fixed, 
the condensates on both sides of the defect will differ only by an Abelian phase. If we 
allow the endpoints to move anywhere on the 5 2 , the phase is 5C/(2)-valued. We study 
the latter in section ||. The configuration is illustrated in figure [l] 

4 Since M = 4 SYM is a conformal field theory, there is also a spontaneous breaking of conformal 
invariance. 

5 This is for a single D(p + 2). For Nf D(p + 2) branes the symmetry is U(Nf). 

6 Although the Dp/D(p + 2) brane configurations with electric flux we have described are supersymmet- 
ric, when they are introduced in the backreacted geometry of the Dp branes they can potentially break 
supersymmetry and have an instability if they are dual to a theory with non-zero chemical potential. We 
want to thank Andreas Karch for pointing out this to us. 
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Dp 



D(p+2) 



Figure 1: Holographic Josephson junction. A D(p + 2) brane is localized at z = and wrapping a 
S 2 . Dp branes end on the D(p + 2), that is seen as a codimension one defect in their worldvolume. 
The separation in the directions of the S 2 between the two Dp halves corresponds to the jump in 
the phase of the condensate through the junction. 

This gives a very simple and geometrical picture of a Josephson junction. Furthermore, 
the degeneracy of ground states makes this a good scenario for studying non-Abelian effects. 
However, because the theory is super symmetric, the system in this form does not admit a 
Josephson effect , since there is no force between the endpoints of the Dp-branes. In order 
to observe a Josephson effect one needs to introduce fluxes on the D(p + 2) brane, as we 
will explain in the next section. In section || we discuss other possible configurations that 
may show a Josephson effect even in the absence of fluxes. 

3. An AC Josephson effect from D-branes 

We now show how to induce an AC Josephson effect in the Dp/D(p + 2) system. We will 
use the following coordinates for the (8 — j?)-sphere 



The D(p + 2) brane wraps the coordinates (j> and 9, while the endpoints of the Dp branes 
are localized in all directions on the (8 — p)-sphere. The directions along the D(p + 2) 
brane which are transverse to the Dp branes are ((f), 9, r). One can then see the endpoints 
of the Dp branes as magnetic monopoles in a 3+1 dimensional spacetime. The Dp branes 
are flat along the x , ■ ■ ■ spatial directions but have a non-trivial profile in time and 

along the x p direction. Thus the only relevant coordinates that span the worldvolume of 
the Dp are x° = Rt and x p = Rz. We will rescale the remaining spatial coordinates and 
the radial coordinate r = Ru by a factor of R in such a way that they are dimensionless. 
The profile for each half- Dp is given by 



dtt\_ p = d9 2 + cos 2 9dn 2 ^ p + sin 2 9d<p 2 . 
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and the rest are zero or trivial. 

The action for each half Dp brane is then 

S Dp = -V P J d p+1 x (e-V=S- Cp+i) 

= -T p R p+1 u 7 - p J d p+1 x (yj 1 + u{f 5 {d^d v <P + dy.ddud + UQ 2 df,5ud v 5u) - 1 

~ -T p R p+1 ul J cF+^rT {d^dA + d^O + u^ 2 d^ud u 5u) + • • • , (3.3) 

where the Dp brane tension is 

Note that the approximated action does not depend on the dimensionality of the dual theory 
except in the overall constant factor. This should not be surprising, as the interpretation 
of (j) a is that they are the Goldstone bosons associated with the U{\) subgroup of the 
spontaneously broken global symmetry of the dual field theory: the free action for such 
fields must take this form. The equations of motion for the linearized fluctuations are just 
the Klein-Gordon equation in flat spacetime 

□0 = 0, ne = o, nsu = o. (3.5) 

The momentum densities are 

tt£ = -T p RP +1 u 2 d^ip, vr£ = -T p R p+1 u 2 d»e, < = -T p R p+1 d»5u. (3.6) 

For the rest of this section, we assume that the radial fluctuations 5u and the the fluctua- 
tions in 9 vanish. The ends of the Dp-branes only fluctuate about their (antipodal) values 
4> = and <f> = 7r. This is the 'classic' Josephson effect, in which the amplitude of the 
condensate on either side of the junction is equal, and the phase is Abelian. 

Given that 4> a correspond to the phase of the condensate, the superfluid current is the 
conjugate canonical momentum 

JM = T M = SS DP = T R P+1 U 2 ^<j) 



From the perspective of the brane it corresponds to components of its energy-momentum 
tensor. For instance, it® is the density of momentum along the (f> direction. In the absence 
of forces 7Ta is a conserved current, but if there is a force in the <fi direction fy, then the 
momentum density current will have non-vanishing divergence 

= f$, (3.8) 

In particular, for homogeneous configurations 

do*° a = f$ d J° a =f$, (3.9) 
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i.e. the time derivative of the superfluid charge density will be non-zero. 

We can induce such a force by using the fact that, from the perspective of the fields 
living on the D(p + 2) brane, the endpoints are magnetically charged objects. They present 
a monopole anti-monopole system. In a supersymmetric configuration, the force induced 
by the bosonic fields on the D(p + 2) cancels against the force induced by fermionic fields. 
In order to generate a non-zero force along the (f> direction, we must introduce a magnetic 
flux oc Bfi, which corresponds to turning on the F$ u components of the field strength 
on the D(p + 2) brane. In general this requires solving non-linear PDEs. We simplify 
the problem by assuming that the displacement of the Dp branes around its equilibrium 
position is small, which is a good approximation as long as the distance between the 
endpoints is not too large. This requires that the magnetic field on the brane is small and 
has an oscillatory behavior. In principle we could also introduce a source in the r direction 
by turning on the Fg^ components of the field strength on the D{p + 2) brane. From the 
point of view of the D(p+2), the endpoints of the Dp are magnetic monopoles. This should 
be analogous to the problem of a string ending on a D-brane (see for instance [29] ) , except 
that the magnetic and electric fields switch roles. Note that the D(p + 2) worldvolume is 
(p + 3)-dimensional, so the magnetic dual field strength 

F = *r+ 3 F, (3.10) 

is a (p + l)-form, and the magnetic dual potential a p-form. This potential has a natural 
coupling to the intersection between the Dp and the D{p + 2), that is p-dimensional. Now 
we turn on the magnetic potential A on the -D(j>+2)-brane. The action on the Dp endpoints 
at z = should have the form 

f dX Ml dX M ' 2 8X m p 

Sa = i„„, dPxA ' hl " M ' {X) — -&r ■ ■ ■ • (3 n) 

with X M the coordinates pullback of the scalar fields on the Dp to the Dp-D(p + 2) 
intersection. We implicitly hide whatever factors of uq etc that arise into the normalization 
of Am\-m v - Since we want a uniform charge density in the x 1 to x p_1 directions, we only 
have nonzero components Amu---p-i where M can be t,u,8,4>. 
Now let us vary the total action 

S (S Dp + S A ) =T P J <P +1 x^5X M UX M (3.12) 

r M \ N ax M > dx M - dx M P - 

+ * dFxbX -T p ^-gg MN g d z X + J r MM 1 M 2 -M p —Q t t^- 



dxP- 

Thus we see that the force on the Dp endpoint acts as a modification of the z = boundary 
condition so that 

. 22 jy dx M ^dx M * dx M v 

T P V-ggMNg o z x = T M M x M 2 -M v — t — ^ QxV _ x ■ (3.13) 



J 2=0 



7 Note that this magnetic field is a flux in the gravity dual and is not a magnetic field in the superfluid 
itself. 



- 8 - 



For instance, if we turn on a magnetic potential of the form A$\2— p(X) = B cos(QX°) 
we would have 

ulT p R p+1 d z (/)(z = 0) = Bsm(Qt) . (3.14) 

and d z 9 = 0, d z u = 0. Electric fields would involve a coupling of the boundary conditions 
for the coordinates <f> and 9, for example. The magnetic field is B = eBq, where e = ±1 
depending on the orientation of the endpoint. 

There is an obvious solution to the equations of motion (|3.5| ) with these boundary 
conditions: 

7T B 

u(t,z) = u , 6{t,z) = -, <f>(t, z) = (j)Q- u<2ji Rp+ln sin(0(t - z)). (3.15) 

This describes a right-moving wave on both Dp branes. A left-moving solution is also 
possible. The charge density transmitted through the junction is 

Q = 9 ttI\ z=0 = -5fisin(fit). (3.16) 

If we add the change in both branes we see that the total variation vanishes. We can then 
compute the differential current density through the junction using charge conservation 

d T z 

— =Bnsm(ttt), (3.17) 
dz 

If we assign a width i to the junction 8 , then the current per unit area across the junction 
would be 

J z = -Q£ = J3JMsin(fit). (3.18) 

This corresponds to an AC Josephson effect with amplitude I c = \B\Q£ and a voltage 
across the junction U = 

We can also describe an approximate DC Josephson effect. First we introduce a mag- 
netic field of the form 

B 4> = Bean(Qt), (3.19) 

so that the charge transfer is 

Q = -BtlcaB(SH). (3.20) 

Now we take the limit of small frequency and large amplitude — >• and \B\ — > oo with 
BO, = I c fixed. Then, the charge transfer becomes approximately constant in time 

Q = -I c [l + 0((nt) 2 )] . (3.21) 

This would be a good approximation as long as tit <C 1. 



8 We introduce this length for the purposes of aligning ourselves with the condensed matter literature. 
In the framework of the model, the junction has zero width - or at most a width near the string scale. 
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4. Non-Abelian Josephson junction 

So far we have studied a case where the difference between the condensates at both sides 
of the junction is just a phase, corresponding to the separation of the endpoints of the Dp 
branes along a circle inside the S 2 wrapped by the D(p + 2) branes. In order to study 
non-Abelian effects we now will allow the endpoints of the Dp brane to move on the full 
S 2 : 

X° = t X* = <f>(t,z) X = 9(t,z) X u = u + 5u(t,z), (4.1) 

where Su will be taken to be small, but <j) and 9 can have variations of order one. We 
will work in an "adiabatic" approximation, meaning that we assume the derivatives of the 
fields to be very small, though some of the field themselves can take on finite values. Then, 
to leading order, the DBI action on the Dp worldvolume becomes 

S Dp = -T P R P+1 J d p+1 x- [ul (sin 2 Od^cj) + d f ,9d fl 9) + d^ud^Su] . (4.2) 

The equation of motion and boundary conditions of 5u are unaffected in the new expansion, 
while the equations of motion of eft and 9 are modified to 

d M [sin 2 9d^(f\ =0, DO- sin 9 cos 9d^(p = 0. (4.3) 
These equations admit solutions that are a superposition of plane waves: 

0(M) = j l^-^M, <Kt,z) = f ^ e M*-*)0( w ). (4.4) 
Therefore, one can introduce boundary conditions with arbitrary time dependence and 

d z 9(t, z = 0) = -d t 9(t, z = 0), d z (j)(t, z = 0) = -d t cj)(t, z = 0). (4.5) 

The boundary conditions follow from the same analysis ( 3. 13| ) . Keeping only the leading 
terms in the derivative expansion we find: 

u 2 T p R p+1 sin 9d z (j) — J 7 (pti2---p— l + J~<f>6i2---p~ idt&i 

u 2 T p R p+1 d z 8 = F et i2- P -i + ^i2- P -M- (4.6) 

The first terms on the r.h.s. of the equations are the magnetic duals to magnetic fields on 
the D(p + 2) F u g and F ur p, while the last terms are dual to an electric field F tu . We will 
introduce the following fluxes on the D(p + 2): 

F(j)ti2-p-i = £?0sin(fit), Foti2-p-i = -Becos(Ot), T<$v2- v -\ = — Fe<p\2- P -i = E. (4.7) 

In general E depends on the u coordinate, but the pullback to the endpoint is a constant. 
We assume that the frequency and the magnetic fields are all small, while the electric field 
is of order one. In terms of a small parameter e, 

d z ~ dt ~ ~ B e ~ f2 ~ e, E ~ 1. (4.8) 
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This way all terms in the equations are of the same order. When ( f4.5|) is satisfied, the 
boundary conditions become 



ulT p R p+1 sin 2 9d t (t) = B^ sin(fti) + Ed t 0, 
ulT p R p+1 d t 9 = B e cos{Qt) - Ed t (ft. 



For simplicity we define rescaled time and fluxes 

E 



so the equations become 



u 2 T p RP+i 



2itB e 
■irr pp+io' 6 



2i:Ba 



u I T p RP+ 1 a 



ulT p RP+ 1 n 1 



(4.9) 



(4.10) 



sin 2 8d T (j) = sin(27rr) + ed T 9, 
d T 9 = bo cos(27tt) — ed T cft. 



(4.11) 



We can find a relation between the change in (ft and in 9 by integrating the last equation 
over a period 

A9= [ drd T 9 = e I dr d T <ft = eAcft. (4.12) 
Jo Jo 

Differentiating with respect to r there are similar relations for higher derivatives, in general 

A(d?0) = eA(d», n = 0,l,---. (4.13) 
Solving for dt<p in the second equation, we get 

8 T (ft = - (8 T 9 + b e cos(2vrr)) . (4.14) 
e 

Plugging this expression in the first equation we are left with the first order equation: 

(e 2 + sin 2 9)d T 9 + b e cos(2vtt) sin 2 9 + eb^ sin(2vrr) = 0. (4.15) 



The integration of ( 4.14 ) leads to 



e \ 2tt J 



(4.16) 



We cannot solve the equations analytically for general e, 6^,6^,, and resort to a numerical 
solutions. An example of a trajectory on the S 2 is given in figure 

When either of the magnetic fields vanishes, we can solve the equation ( [4. 15 ) analyti- 
cally (up to a transcendental relation). When be = 0, we have 



2 1 
e 2 + - 



- sin(20) = cos(2vrr) + C 

4 2ir 



e 



+ 



1 



- sin(20) - %± COS {2tt9) + C 
4 2tt 



(4.17) 
(4.18) 
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""6.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

<f>/7T 

Figure 2: The line represents the evolution of the solution in the (</> / V , 9 / V) plane, starting at 
(1,1/2). After one period where the magnetic fields return to the same values, the value of the 
scalars is different. The values of the parameters for this solution are e = —1, be = —5 and = 5. 

where Co is the initial value of the transcendental, 

C = (e 2 + 0(r = 0) - ~ sin(20(r = 0)) + ^. (4.19) 

The transcendental equation has only one solution for each value of r. As one can see from 
the periodicity of the solutions, then, in the case of bg = 0, both eft and 6 return to their 
original values after one cycle. The same is true for b^ = 0, where the solutions for the 
phase are 

(9 -e 2 cos (9 = -— sin(2vrr) - D (4.20) 
2tt 

(t> = <t>o + -(D + e 2 cot9) (4.21) 

where 

D = 6{t = 0) - e 2 cos(#(r = 0)) . (4.22) 

We can now see that in order to measure a non-trivial Berry holonomy in the system, we 
must turn on b^, be, and e. 

5. Berry holonomy and connection 

In the holographic Josephson junction there is a degenerate space of configurations, that 
corresponds to moving the endpoints of the Dp branes to different positions on the S 2 
wrapped by the D(p + 2) brane. We have seen that as a consequence there is a non- 
Abelian Josephson effect and that the final state of an adiabatic evolution of the system 
along a closed curve in the space of magnetic fields is different from the initial state. We 
can make this statement more quantitative by defining a Berry connection in the space of 



l.U 



0.8 



0.( 



7T 



0.4 



0.2 



n ii 
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magnetic fields and measuring the Berry holonomy along the closed path. As we will see, 
the evolution of the phase of the condensate and the Josephson current is determined by 
parallel transport along the curve with respect to the Berry connection. 

Let us start mapping the two-sphere to the complex plane using a stereographic pro- 
jection 

z = tan^e i *. (5.1) 

The group of automorphisms of the sphere is GL(2, C), i.e. 2x2 complex matrices with 
non-zero determinant: 

'-(:;)■ <»> 

The action of the group on the complex plane is 

' " aZ + b /r o\ 

z =gz = — — . (5.3) 

cz + a 

We will actually be interested in a SU(2) subgroup. 

In the model of non-Abelian holographic Josephson junction we are studying, the mag- 
netic fields are external parameters that we vary adiabatically. As the magnetic fields vary, 
the endpoints of the Dp brane move on the two-sphere. This is seen as a change in the state 
of the system, that we can map to a trajectory in the complex plane through the stereo- 
graphic map. This trajectory can be described using the group of automorphisms above. 
For instance, if the initial and final points of the trajectory in one period of oscillation are 

Zj=tan^, z f = tan (t+l£^ ^ ■ ( 5 - 4 ) 

The two are related by the action of an element of an SU(2) group 
Zf = g(a f ,Pf,0)zi, 

( e ia fl 2 \ / cos sin ^ \ 
9 (a f J f -0)= 9 ,(a f ) 9e m =1^ Q e _° //2 J (^ n | j . (5.5) 

If the initial point is not real Zi = tan |e*^, then the group element relating the initial and 
final points is modified to 

g(a f ,P f ;cl>)=g^a f )g (t> (cl>)ge(Pf)g 4> (-<f>)= ^ e _° //2 j ^ _ j^-/Tos4f )" 

2 (5-6) 

We can divide the trajectory into N small pieces characterized by the time intervals 
At = ti + i — U, i = 0, • • • , N. When they are glued together it is clear that the group 
element that relates the initial and final points is the product of group elements that relate 
the endpoints of each of the smaller intervals: 

N-l 

9 {af,Pf-m) = II gWu),P(u);Hu))- (5-7) 

i=0 



- 13 - 



We can make this an statement about an infinitesimal change along the trajectory. The 
change in the complex plane along the trajectory is 

z(t + At) - z(t) = d t zAt = - ( 1 + tan 2 - ) e^d t 9At + itan -e i(t >d t (f)At. (5.8) 



On the other hand, a transformation ( |5.3[ ) with infinitesimal values of a and j3 in ( |5.6| ) is 

1 / 9\ 9 

gz - z = -11 + tan 2 - + i tan -A. (5.9) 



Therefore, we can identify (3 = dtOAt = 9 At and a = dt4>At = <pAt. We can describe the 
trajectory using the equation 

z = ih(5a, 5b, 5c, 6d)z, (5.10) 
where h is an element of the gl{2,C) algebra, acting on z as 

h(5a, 5b, 5c, 5d)z = (5a — 5d)z + 5b — 5cz 2 . (5-11) 

In our case it is an element of the su(2) subalgebra. Although we do not know h a priori, 
we can extract its value from the solutions we have found: 

3 2 1 

h{t) = 0y + 0cos0y +^sin0y, (5.12) 

where a 1 are the Pauli matrices. 

In this way we can associate to each point of the trajectory an element of the su(2) 
algebra. We can rewrite ( |5.10| ) as the equation of parallel transport along the trajectory 
in parameter space. This allows us to define an su(2) connection along the closed curve 
in the space of magnetic fields. We will identify this with the Berry connection along the 
curve. The unit vector tangent to the curve is 

' (OtB^dtBg). (5.13) 



y/{dtB^ + {d t B e y 
Then, equation ( |5.10D becomes 

t-(V-iA)z = 0, (5.14) 
where the value of the Berry connection along the curve is then defined as 

t-A= - Ht) (5.15) 

^(d t B^Y + (d t B y 

Note that fl5.14| ) determines the parallel transport, as defined by the Berry connection, 
of the phase of the condensate along the curve in the space of magnetic fields. The Wilson 
loop along the closed curve in the space of magnetic fields determines the transformation 
( |5.G| ) that relates the two endpoints of the trajectory in the space of values of the condensate 

g{a f ,P f ,m)=e l§ c t - A . (5.16) 
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By Stokes' theorem, this should be related to the integral of the Berry curvature in the 
area enclosed by the curve 



In our examples the Wilson loop is non-trivial, indicating that there is a non-zero Berry 
curvature. 

6. Conclusions and future directions 

In this paper we have used D-branes to construct a simple holographic model for Josephson 
junctions in p + 1-dimensional superfluids, < p < 4, as a Dp/D(p + 2) intersection. 
By varying the magnetic fields on the D{p + 2) brane we have found a non-trivial Berry 
holonomy, which measures the amount of electric flux, or equivalently the number of strings 
on the D(p + 2) brane. This should be an integer number, so the Berry holonomy must 
be quantized. It may thus be used as a way to characterize new topological phases of 
holographic superfluids. 

For p = 3 one could define a different kind of topological phase if one restricts to a 
particular class of models. The starting point is a single D5. By adding a large number of 
A^S'5-branes one can produce a super symmetric Janus geometry where the theta angle on 
the D3 jumps across the defect, as suggested in [23]. In this case a Chern-Simons term for 
the D3 brane gauge fields is induced on the intersection between the D3 and the 5-branes. 
This implies that there is a Hall effect on the defect, which one can see as the boundary 
between two systems. 9 This is quite similar to the effective theory that describes the re- 
sponse of topological insulators to external electromagnetic fields [26,27]. In this sense one 
can see the coefficient of the Chern-Simons term as labelling distinct "topological phases" , 
its value is quantized since the numbers of D5 and NS5 branes should be integers. 10 In the 
absence of NS5 branes there is no Chern-Simons term, so this would correspond to the 
"trivial phase" . When the fluxes on the D5 brane are turned on, the boundary conditions 
of the fields on the D3 brane change. The D3/D5 intersection with string flux also falls into 
the supersymmetric cases studied in [22], so it would be interesting to see the connection 
of the boundary conditions with the Berry holonomy. 

There are many other interesting directions which remain to be explored. 

An obvious extension is to study intersections with more branes, n p Dp and n p+ 2 
D{p + 2) branes. The relevant group of symmetries at the intersection is enlarged to 
U(n p ) x C/(n p +2) x SO(3), which can be broken in different ways as the Dp branes are 
split, leading to different types of Josephson junctions. With several D(p + 2) branes it 

9 A Hall effect can also be induced if there is a background axion and magnetic flux on the sphere that 
the D5 is wrapping [25]. 

10 This theory is not unique. In principle one can construct many supersymmetric Janus geometries where 
the theta angle is not quantized, but those would not correspond to a distribution of NS5 branes. We thank 
Andreas Karch for illuminating comments on this issue. 




(5.17) 
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is also possible to build arrays of Josephson junctions, by separating the branes along the 
transverse spatial direction and connecting them with Dp branes. 

Another clear extension of this work would be to study other Josephson junctions 
obtained from .D-brane configurations. A naive candidate are p-dimensional intersections of 
Dp and D(p+4) branes. However, although in this configuration there is a codimension one 
defect on the Dp branes, they cannot be split in two halves ending on the D(p + 4) branes. 
For p = 3 there are no other possibilities left. For p = 1,2 the remaining possibilities are 
Dp/D(p + 6) intersections D1/D7 and D2/D8, which are also super symmetric. 

Although the simplest models are supersymmetric, slightly more complicated models 
could lead to non-supersymmetric Josephson junctions, with a non-zero Josephson current 
even in the absence of additional fluxes on the defect. A possible candidate is the (2 + 1)- 
dimensional intersection D3/D7' studied in [30] as a holographic model with a Quantum 
Hall Effect. Although the D3 cannot end on a D7, the D7' has fluxes turned on that can 
be seen as D5 branes dissolved in its worldvolume, on which the D3 brane can in principle 
end. 

Supersymmetry could also be broken explicitly, for instance by making one of the 
spatial directions along the D(p + 2) compact with length L and imposing antiperiodic 
boundary conditions for the gaugino on the D(p + 2): 



If we choose the length of the compact direction to be much smaller than the radius of the 
S 8 ~ p , we can neglect the effect of the massive gaugino on the D(p + 2). L has still to be 
much larger than the string scale in order to stay in the supergravity approximation. This 
does not affect the stability of the D(p + 2) brane at the classical level, which is dependent 
only on the bosonic fields. The final picture is that we put in contact two superconductors 
of the same material but with different condensates along a strip of width L. 

One could also break supersymmetry by modifying the background geometry, for in- 
stance by introducing a temperature. This would also be interesting in this framework in 
order to be able to plot the phase diagram. However, if one naively replaces the Dp-brane 
background with a black hole, the probe branes simply fall in, so the solution is not stable. 
One would therefore require a more complicated background, including additional bulk 
fluxes to counter the effect of gravity. Most likely in order to be able to pull Dp branes 
outside the horizon it would be necessary for the background itself to be dual to a super- 
fluid phase, so the phase diagram will be determined by the background. There are several 
models constructed from consistent truncations of supergravity (e.g. [33-37]), although in 
order to introduce D-branes one should uplift them to ten dimensions. 

Supergravity backgrounds can also be useful to go beyond the probe approximation. 
An approach directly related to the D-brane constructions we have presented will be to 
replace the defect brane by a Janus geometry [31,32], while keeping the Dp branes as 
probes. For the intersection of D3 branes with 5-branes there are known solutions [23,38] 
which include configurations with the 5-branes separated in different stacks [39,40]. 





(6.1) 
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